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We determine the class of damped modes y which are related to the common free damping modes y by supersymmetry. 
They are obtained by employing the factorization of Newton's differential equation of motion for the free damped oscillator 
by means of the general solution of the corresponding Riccati equation together with Witten's method of constructing the 
supersymmetric partner operator. This procedure leads to one-parameter families of (transient) modes for each of the three 
types of free damping, corresponding to a particular type of antirestoring acceleration (adding up to the usual Hooke restoring 
acceleration) of the form a{t) — ^^^^j^^a y, where 7 is the family parameter that has been chosen as the inverse of the Riccati 
integration constant. In supersymmetric terms, they represent all those one Riccati parameter damping modes having the same 
Newtonian free damping partner mode. 
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The damped oscillator (DO) is a cornerstone of physics 
and a primary textbook example in classical mechanics. 
Schemes of analogies allow its extension to many areas of 
physics where the same basic concepts occur with merely 
a change in the meaning of the symbols. Apparently, 
there might hardly be anything new to say about such an 
obvious case. However, in the following we would like to 
exhibit a different and nice feature of damping resulting 
from the mathematical procedure of factorization of its 
differential equation. In the past, the factorization of the 
DO differential equation (Newton's law) has been tackled 
by a few authors [Q but not in the framework that will 
be presented herein. Namely, recalling that such factor- 
izations are common tools in Witten's supersymmetric 
quantum mechanics |^ and imply particular solutions of 
Riccati equations known as superpotentials, we would 
like to explore here the factoring of the DO equation by 
means of the general solution of the Riccati equation, a 
procedure that has been used in physics by Mielnik ||] 
for the quantum harmonic oscillator. In other words, our 
goal here is to exploit the nonuniqueness of the factoriza- 
tion of second-order differential operators, on the exam- 
ple of the classical damped oscillator. By doing this one 
may hope to gain insight into the free damping motion. 
We write the ordinary DO Newton's law in the form 



differential equation 



Ny^[^ + 2(3^^+f3-]y=ip--cu-)y 



a^y 



(1) 



Ly= \ P\y± = ±ay± 



(3) 



It follows that y± = g-l3t±at g^j^j ^^^g ^g^^^ build 
through their superposition the general solution as 
y = e^^*{Ae"* + i3e~"*). Thus, for free underdamp- 
ing, the general solution can be written as j/„ — 
Ae'l^'^ cos{\/-a'^t + (p), where A = 2^/\AB\ and (j> = 
Arcos(^^i^), whereas the overdamped general solution 

is Ae~'^^cosh{at + (p), where A = 2^\AB\ and <j) = 
Arcosh(^^i^). The critical case is special but well known 
having the general solution of the type yc ~ e^^*{A + 
Bt) . 

Since, as we mentioned, the factorization given by Eq. 
(2) may not be the only one possible, let us now write 
the more general factorization 
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where f{t) and g{t) are two functions of time. The con- 
dition that Ng be identical to N leads to f{t) -\- g{t) = 2/3 
and g + fg = (3, that can be combined in the following 
Riccati equation 



-f -f + 213 f = (3^ 



(5) 



i.e., we already added a 0^y term in both sides in order to 
perform the factoring. The coefficient 2(3 is the friction 
constant per unit mass and is the natural frequency 
of the oscillator. The factorization 



a^y 



(2) 



follows, and previous authors |^ discussed the classical 
cases of underdamping [a^ < 0) , critical damping (a^ = 
0), and overdamping (a^ > 0) in terms of the first-order 



By inspection, one can easily see that a first solution 
to this equation is f{t) = f3 {g{t) — (3), which is the 
common case discussed by all the previous authors [0. 
Changing the dependent variable to h{t) = f{t) — (3, we 
get a simpler form of the Riccati equation, i.e., h [t) + 
— 0, with the particular solution h{t) = 0. However, 
the general solution is h{t) = 



as one can 



t+T 7*4-1 ' 

easily check. The constant of integration T = I/7 occurs 
as a new time scale in the problem; see below. Therefore, 
there is the more general factorization of the DO equation 



I 



than Eq. (2) 
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A few remarks are in order. While the hnear operator 
L = ^ + /3 has y± as eigenfunctions with eigenvalues ±a, 
the quadratic operator N has y± as degenerate eigen- 
functions, with the same eigenvalue a^. On the other 
hand, the new linear operators and do not have 
y± as eigenfunctions since A'^y± = {±a + 7^^jj)y± and -^(j-^ 



A y± ~ (±a — ^^^i )y±, although the quadratic oper- 
ator Ng = A'^A~ still has y± as degenerate eigenfunc- 
tions at eigenvalue a^. We now construct, according to 
the ideas of supersymmetric quantum mechanics ||] , the 
supersymmetric partner of Ng 



Ng = A-A-^ 



27^ 



(7i + l)2 



(7) 



This second-order damping operator contains the addi- 
tional last term with respect to its initial partner, which, 
roughly speaking, is the Darboux transform term ^] 
of the quadratic operator. The important property of 
this operator is the following. If y^ is an eigenfunc- 
tion of Ng^ then A^j/o is an eigenfunction of Ng since 
NgA-yn A-A+A-yo = A- Ngyo and Ngyo = a^yo, 
implying Ng{A-yo) = A- Ngyo = a'^iA-yo). The con- 
clusion is that Ng has the same type of "spectrum" as 
Ng, and therefore as N. The eigenfunctions y± can be 
constructed if one knows the eigenfunctions y± as 



y± = A y± 
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di 
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and thus 



y± 



7< + 1 



-I3t±at 



(8) 



(9) 



These modes make up a one-parameter family of damping 
eigenfunctions that wc interpret as follows. We write 
down the usual form of the Newton law corresponding to 
the Newton operator Ng, 



..d 



27^ 



(7< + l)2 



y = o 



(10) 



Examination of this law shows that the term 27^/ (7t -t- 
can be interpreted as a time-dependent antirestoring 
acceleration (because of the minus sign in front of it) 
producing in the transient period t < 1/(3 the damping 
modes given by y above. 

We present now separately the y families of modes cal- 
culated as superpositions of the modes y± for the three 
types of free damping. 

(i) For underdamping 0^ < let a = itoi, where 
uji — y^ujQ — l3^. The original eigenfunction is yu = 



A.a cos{uJit + 4))e while the supersymmetric family is 
iju = -Au[u}i sin(wit + 4>) + cos(a;it + 4>)]e~'^K 

(ii) In the case of critical damping = Wg, the general 
free solution is yc = Ae~f^* + Bte~^* , whereas the tilde 
solution wiU be y^ = + ^{-ft + l)^]e-'^*. There is 

a difficulty in this case since y+ = A^y^ = "^^^^ e"^*, 
whereas y_ — A^y^ — ^^^-^ e~^* cx y+. To find the 
independent y_ solution we write y^ = z(i)y+ and de- 



termine the function z{t) from Ngij^ ~ 0. The result is 
where C is an arbitrary constant, and 
£) (7*-Hi) g"/3t^ £) being another arbitrary 



therefore y_ 
constant. 

(iii) For overdamping 0^ > Wq, the initial free general 



solution is yo — AoC cosh(at -|- 
tion is yo — ~Aoe~^^[asmh.{at + 



whereas the 7 solu- 
-:^cosh(at-|-0)]. 
Plots corresponding to these cases are presented in 
Figs. 1-3. We note that in the limit 7^0 the modes 
yuj,i3,-y are going to the Newtonian damping modes y^^^^ 
for all three classes of free damping motion. Moreover, we 
placed ourselves herein in the well-behaved regime of mo- 
tion, i.e., for time and parameter ranges where the modes 
do not grow with time and their amplitudes are finite. 
However, from the point of view of the 7 parameter the 
modes y are always singular, i.e., they blow up at some 
negative time moment for positive 7 and at some positive 
instant for negative 7. Such blow-up solutions are quite 
well known in nonlinear physics. On the other hand, 
even the Newtonian modes yu},[3 are growing with time 
in the past or for negative /3 in the future (divergent and 
flutter instabilities are textbook knowledge [D). What 
we claim here is that when one starts a damping-type 
measurement after a "mechanical" blow-up phenomenon, 
Riccati parameter modes may be present. As we said, 
they may also occur before a blow-up phenomenon (for 
negative 7), an equally important case. In this situation 
the Riccati parameter distinguishes them from more com- 
mon instability modes. Thus, an extended, Riccati-type 
parametrization of free damping can indeed be useful. 
The complexification of the Riccati parameter adds one 
more parameter to the Riccati damping modes. Depend- 
ing on the sign of the imaginary part, new contributions 
to either damping or destabilization of the modes occur. 

In summary, what we have obtained here are Riccati 
parameter families of damping modes related to the New- 
tonian free damping ones by means of Witten's super- 
symmetric scheme and the general Riccati solution. 
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FIG. 1. Initial 

free underdampcd mode of the type yu ~ e~*^^''cost (bold 
curve) and members of its 7 family of supcrsymmetric damp- 
ing modes yu = — e~'''^°(sinf -|- ^^^^cost) for the following 
values of parameter 7: dashed curve - 1; bold dashed curve - 
1/2; solid curve - 1/10. 
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FIG. 2. Initial free critical damping mode yc = e '(1 + 1) 
(bold curve) and members of the corresponding 7 family 

yc = + ^''^^i^ ) for the 7 parameter taking the 

following values: dashed curve - 5; bold-dashed curve - 
5/3;dot-dashed curve - 1. 



FIG. 3. Initial free overdamped mode of the type 

j/o ~ e~* cosh(t/5) and members of its supersymmetric 7 fam- 
ily §0 = e *[-| sinh(t/5) — cosh(t/5)] for the following val- 
ues of the parameter 7: dashed curve - 1; bold dashed curve 
- 1/2; solid curve - 1/10. 
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FIG. 4. The antirestoring acceleration for the under- 
damped modes at the values of the 7 parameter: (a) 1; 
(b) 1/2; (c) 1/10. 



FIG. 5. The antirestoring acceleration for the critical 
modes yc at the values of the 7 parameter: (a) 5; (b) 5/3; 
(c) 1. 
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FIG. 6. The antirestoring acceleration for the ovcrdamped 
modes yo at the values of the 7 parameter: (a) 1; (b) 1/2; (c) 
1/10. 
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